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Abstract. We derive a new Banks-Casher-type relation which relates the density of complex Dirac eigen- 
values at the origin to the BCS gap of quarks at high density. Our relation is applicable to QCD and 
QCD-like theories without a sign problem, such as two-color QCD and adjoint QCD with baryon chemical 
potential, and QCD with isospin chemical potential. It provides us with a method to measure the BCS 
gap through the Dirac spectrum on the lattice. 
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1 Introduction 

Symmetries and their spontaneous breaking are impor- 
tant ingredients for a modern understanding of physics. 
In quantum chromodynamics (QCD), the spontaneous 
breaking of chiral symmetry, characterized by the chiral 
condensate (ipip) > gives rise to a variety of phenomena in 
low-energy hadron physics. Understanding the nature and 
the mechanism of chiral symmetry breaking is a funda- 
mental challenge for theoretical approaches to QCD. 

More than 30 years ago, it was realized by Banks and 
Casher [1] that in the QCD vacuum the Dirac eigenvalue 
spectrum at the origin is related to the magnitude of the 
chiral condensate in the chiral limit through the relation 



(1.1) 



where p(X) is the spectral density of the Dirac operator. 
This relation allows us to measure the magnitude of the 
chiral condensate through the Dirac eig enva lue spectrum 
in lattice simulations [2,3,4 . Relation (1.1 1 was derived 



for purely imaginary (or real, depending on the conven- 
tion) eigenvalues. It is valid not only in the QCD vac- 
uum but also at nonzero temperature T, as T does not 
invalidate any part of the original derivation. Since p(0) 
is nonzero in the phase with broken chiral symmetry, 
whereas it vanishes in the high-temperature chirally sym- 
metric phase, it serves as a clean indicator of the chiral 
phase transition 5 , 6 , 7 J . In this way the Dirac eigenvalue 
spectrum plays a crucial role in the modern understand- 
ing of the nonperturbative phenomena of chiral symmetry 
breaking and restoration. The Banks-Casher relation is 



also valid in QCD-like theories, including two-color QCD 
and QCD with fermions in the adjoint representation, 
which share some universal features with bona fide QCD 



An interesting question to ask is whether the Banks- 
Casher relation could be generalized to QCD and QCD- 
like theories at nonzero quark chemical potential //, whose 
phase diagrams have been studied extensively (see, e.g., 
[mfT2] for recent reviews). In fact, another interesting 
nonperturbative phenomenon is expected to occur at large 
p |13j : quark-quark pairs are formed near the Fermi sur- 
face due to the Bardeen-Cooper-Schrieffer (BCS) mecha- 
nism. Connections between BCS pairing and the Dirac 
eigenvalue spectrum were first investigated in [14 for 
three-color QCD and in [T5l[T6l[T7] for two-color QCD. 
Nevertheless, the Banks-Casher relation itself has not yet 
been generalized to dense quark matter)]] 

In this paper, we derive a new Banks-Casher-type rela- 
tion in theories without a sign problem, which states that 
the BCS gap (or the superfiuid gap) of quarks at asymp- 
totically high densitjf^] in the chiral limit is proportional 
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A Banks-Casher-type relation that relates the spectrum 
of the Dirac singular values at the origin to the condensate 
of quark-quark pairs in QCD-like theories at any fi has been 
obtained in |18| . 

2 Strictly speaking, the new Banks-Casher-type relation is 
not exact at intermediate density. For example, instanton ef- 
fects, which are negligible at asymptotically high density, may 
quantitatively alter our relation at intermediate density (see 
footnote[7|. Still, instanton effects are suppressed by high pow- 
ers of Aqcu/p [T9"ll2"0"] . and are expected to be minor for fj, 
larger than Aqcd- Whether the Dirac spectrum for a given 
H is related to the chiral condensate or the BCS gap can be 
judged by looking at the spacing of the near-zero eigenvalues as 



2 
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to the average density of near-zero complex eigenvalues 
of the non-Hermitian Dirac operator. The explicit form of 
the relation is 

2?r 3 

-Pi(0), (1.2) 



A 2 



rep 



where A is the energy gap in the spectrum of the 
quasiquarks above the ground st ate, pi (A) is the two- 
dimensional spectral density, see (2.8) for the definition, 



and G? rcp is the dimension of the color representation in 
which the quarks transform: c?f un d = N c and d a( jj = AT? — 1 
for the fundamental and adjoint representation of SU(iV c ), 
respectively. We prove this relation for two-color QCD 
[2THTU] with baryon chemical potential, QCD with isospin 
chemical potential [2211231124] . and QCD with fermions in 
the adjoint representation |10j with baryon chemical po- 
tential All of these theories have a positive-definite path 
integral measure if the parameters of the theory are chosen 
properly. On the other hand, our relation does not extend 
to color-superconducting phases in three-color QCD with 
baryon chemical potential]^] because the spectral density 
of the Dirac operator near the origin seems to be well- 
defined in the thermodynamic limit only for theories with 
a positive-definite path integral measure [8, 27, 28 (see also 
appendix C of [T5] for a detailed discussion of this point). 
Nevertheless, our new relation serves as a useful means of 
measuring the magnitude of the BCS gap at high density 
in lattice simulations of theories without a sign problem. 
Together with the rigorous spectral sum rules for Dirac 
eigenvalues in the microscopic limit |15j . our result con- 
stitutes a unique and quantitative way of studying the 
physics of relativistic dense matter. 

This paper is organized as follows. In sections [2j [3j and 
[4] we derive the Banks-Casher-type relation for the BCS 
gap in two-color QCD with baryon chemical potential, 
QCD with isospin chemical potential, and adjoint QCD 
with baryon chemical potential, respectively. Conclusions 
and outlook are presented in section [5j In appendix |A"] we 
give the microscopic derivation of a high-energy constant 
appearing in the effective Lagrangian, and in appendix |B] 
we discuss the consequences of a high-energy term that 
was missing in [15j . In appendix[C]we derive the free Dirac 
spectrum at nonzero ji, while in appendix [D] we comment 
on an alternative derivation of the Banks-Casher-type re- 
lation based on the replica trick. 

In this paper we always work in Euclidean spacetime. 
Our convention is that the average of an observable O 



a function of the four- volume V4: If it is proportional to 1 /V>4 
(I/V4), the Dirac spectrum is related to the BCS gap (chiral 
condensate). A deviation from the l/yVj scaling at interme- 
diate density due to instanton effects is expected to become 
smaller as /_i becomes large compared with Aqco- 

3 The Dirac operators in these theories belong to three dif- 
ferent symmetry classes [2">il'2(.i | . 

4 It is worth stressing that a link between the Dirac spec- 
trum and the BCS pairing does exist even for a color- 
superconducting phase [13] . However, in such a case A 2 re- 
ceives contributions from an extended region within the Dirac 
spectrum and is not simply proportional to pi(0). 



over gauge fields weighted by the pure Yang-Mills ac- 
tion is denoted by (0) Y m> wn il e the average including n 
dynamical flavors is denoted by (0) N _ n . By definition, 

(O)YM = (°>iV / =0- 



2 Two-color QCD 

In this section we consider two-color QCD at high den- 
sity, i.e., for sufficiently large quark chemical potential fi 
such that Asu (2) ^ Mi where ylsu(2) is the typical scale of 
the strong interactions in two-color QCD. In section |2.1| 
we derive a Banks-Casher-type relation between the BCS 
gap squared and the Dirac eigenvalue density at the ori- 
gin, working with two flavors for simplicity. In section |2.2| 
we show that this result is consistent with earlier results 



[17] on the microscopic Dirac spectrum. In section 2.3 we 
derive consistency relations based on known spectral sum 
rules 115 . In section 2.4 we show that our results are also 



valid for general even Nf 



2.1 Derivation of Banks-Casher-type relation 



The partition 

by 



function of two-color QCD at \i 7^ is given 



Z({m f }) 



N f 

n det (-° • 



/=i 



(2.1) 



YM 



where D = D(fi) — ^ V D V + /J74 is the Dirac operator and 
the m/ are the quark masses. In the following we take 
Nj — 2 for si mplic ity (comments on general Nf will be 
given in section 2.4 1. The partition can then be written as 



Z(mi,m 2 ) — (det(D + mi) dct(Z? + m 2 )) YM 
= (det(D + mi) det(D f + m 2 )) YM 



(2.2) 



The last line follows from the anti-unitary symmetry 
[Cj^K, D] = of the two-color Dirac operator [5J (with 
the charge conjugation operator C, the generator t 2 of 
the color group, and the complex conjugation operator 
K), which implies D* = (C^t^D {Cj 5 t 2 ). 

To derive Banks-Casher-type relations it is mandatory 
that the integration measure is positive definite, see ap- 
pendix C of [18j . Looking at (2.2), we have two options: 
777,1 = z and 7772 = z* with z complex, or mi = m 2 = m 
with to real. In the following we choose the first option 
and comment on the second one a t the end of this section. 
Before proceeding we note that (2.2) with mi = z and 
7772 = z* is similar in spirit (although not quite identical) 
to the "Hermitization" of a non-Hermitian problem, see 
[2^[3ll3Tll3"2"l[53l^4] for other applications of this method. 

Let us consider the quark mass dependence of the par- 
tition function (or the shift of the free energy due to the 
quark mass) in two-color QCD at high density. Recall that 
there is no contribution of 0{m 1 ) at high density where 
the chiral condensate is vanishing. The lowest-order shift 
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in the free energy is thus 0(m 2 ). In [15] the contribu- 
tion S£ of the BCS pairing of quarks near the Fermi sur- 
face to the mass dependence of the ground-state energy 
density was computed based on high-density effective the- 
ory (HDET) with quark masses [35] . However, this result 
needs a small modification to allow for a mass matrix of 
the form M = diag(z,z*). To see this we first note that 
usually one introduces complex quark masses by writing 
the mass term as Mlr 1 t ! 75 + Mm 1 ~ 75 and then letting 



M LR = M and M RL = Aft. However, in rtO) we have 
chosen Mlr — Mrl = M = diag(roi, TO2), where mi and 
m 2 could be complex^] In this case the result of 15] for ddlog Z(z, z*;e) = ( J]] 
Nf = 2 becomes ' ' * 



introduce a small regulator e > (see [29,30,32,39,40] for 
a similar trick) to obtain 

Z( Zl z*;e) = (l[[(\ n + z)(\* n + z*)+e 2 ]) . (2.6) 

\ „ / YM 

The necessity of this regulator is not obvious at this point 
and will be explained after (2.13). Here we only note that 
it ensures that the factors in the product can never become 
zero. From (2.6 1 we derive 



(|A n +*| 2 + £ 2)5 



8£ — man 

A 



2tt 2 



A 2 (A 2 + A* 2 )detM 



(2.3) 



where A £ U(l). Note that for 5£ to be real, detM has 
to be real, which is the case for M — diag(z, z* ). For this 
choice of M the RHS is minimized for A = ±1. In 0(m 2 ) 
there is also a high-energy ternfj H 2 ^(MrlMlr) = 
H 2 tr(M 2 ) that was not consideredin [T5] so thaiQ 



logZ{z,z*) = — A 2 zz*+H 2 (z 2 + z* 2 )+0(\z\ 3 ). 



where V4 is the four- volume and H 2 is a high-energ y 
stant whose origin will be clarified in appendix AT 
noting d = d/dz and d = d/dz* we obtain 



(2.4) 

r con- 
De- 



lim lim 

z— >o v 4 — >oo 



dd log Z(z, z* 



;A 2 



(2.5) 



i.e., we can extract A 2 by suitable derivatives with respect 
to the complex quark masses. 

We now make contact with the Dirac spectrum by writ- 
ing Z(z,z*) in terms of the eigenvalues A„ of the Dirac 
operator]^] Starting from (2.2) we follow [5S] eq. (56)] and 



5 See also footnote 2 of [36] for a discussion of this subtlety. 

6 We follow the convention in [37]. Our high-energy term is 
tr(M 2 ) instead of the more conventional tr(M^M) for the same 
reason as discussed before ( |2.3[ ). Note that tr(M 2 ) is real for 
M = diag(z,z*). 

7 The explicit breaking of the Ua(1) symmetry due to instan- 
tons is ignored here, based on the fact that at asymptotically 
high density the QCD running coupling is small and instantons 
are screened. In other words, Ua(1) is assumed to be sponta- 
neously broken by (tptp) 7^ 0. If we were to take the explicit 
breaking of Ua(1) into account, yet another mass term of the 
form Ci ns t det M would appear in the free energy, making it dif- 
ficult to extract A 2 separately from Ci nst . Finally, note that the 
contribution from the Bedaque-Schafer term [38] , which affects 
the ground-state energy density at 0(m 4 ), is also ignored here. 

8 Note that the H2 term, which originates from the bulk 
Fermi sea and the vacuum, does not appear in (2.3) because 



HDET captures the physics only near the Fermi surface. This 
is why we have to compute H2 separately from microscopic 
QCD, see appendix [X] 

9 In principle the Dirac operator can have exact zero modes 
with A n = 0, but since at high density topology is suppressed 
we can neglect exact zero modes in this paper. 



E 



E 



A„ + z 



\\ rn + z\ 2 +e 2 ^ |A„ + z| 2 +e 



N,=2 



Nf—2 



\ ^ A n + z* \ / x ^ A„ + z 

^ \X n + z\ 2 +s 2 / Nf=2 \^ \Xn + z\ 2 + e 2 / Nf=2 ■ 

(2.7) 

We now take the thermodynamic limit (V4 — > 00) and 
then the chiral limit (z — > 0). Introducing the one-point 
function^ 



MA) = Jim lim ^(£* 2 (A-A 



z-S-0 VWoo V4 

and the connected two-point function 
p 2 (\,X) = lim lim 

z— yO V4— >oo 



N t =2 



^-(^5 2 (A-A m )^ ( 5 2 (A'-A„)) Ar=2 -y 4 p 1 (A)p 1 (A') 



V A 



we obtain 



1 - f £ 2 
lim lim — 99 log Z(z, z*; e) — / d 2 \-r—z ^r— 

z^OV^ooVi X ' J (|A| 2 +£ 2 ) 2 



(2.9) 



Pi(A) 



d 2 \ / d 2 X' 



X* 



X' 



|A| 2 + e 2 |A'| 2 +e ; 



p^(A,A'). (2.10) 



From chiral symmetry we have p 2 (X,X') = p 2 (—X,X r ) = 
p 2 (X.—X') = p 2 (—X,—X'), and therefore the second line 
vanishes. Using the formula 



5 2 (z) = 1 lim 



(2.11) 



7T a^O (\z\ 2 + a 2 ) 2 

for the delta function in the complex plane we arrive at 
lim lim lim — ddlog Z(z, z*; e) = 7rpi(Q) . (2.12) 

e— >0 z— >0 V^— too V4 



Comparing this to (2.5) we obtain our main result 

a 2 = 4pM, 



(2.13) 



Note two small changes in conventions compared to [15] . 
There, the Dirac eigenvalues were defined to be iX n , and the 
one-point function was defined without the factor of I/V4. 
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which is a Banks-Casher-type relation for the BCS gap 
squared in terms of the Dirac eigenvalue density at zero. 

Our derivation assumes the high-density BCS-type su- 
perfluid phase of QCD. At low density, the Dirac spectrum 
has also been analyzed by means of low-energy effective 
theories at the mean field level [4"Tll4"2"l[2"5] . These analy- 
ses at low density, which complement our analysis at high 
density, showed that the Dirac spectrum at small fi forms 
a band of width 4F 2 fj, 2 /S and height S 2 / (Ann 2 F 2 ) , which 
extends infinitely along the imaginary axis rj Here F and 
S are the usual low-energy constants (pion decay constant 
and chiral condensate) at /i = 0. 

In the rest of this subsection we make a number of 
technical comments that can be skipped in a first reading. 

First, one can show that the contribution of the large 



eigenvalues to the first integral on the RHS of (2.101 is 
suppressed if we introduce an ultraviolet cutoff A\ 
take the limit e — > before the limit /Itjv — > °o- 



itjv and 
Also, 

even if we had not neglected exact zero modes, their con- 
tribution to the integral would have been suppressed in 
the limit V4 — ► 00 in analogy to the discussion in [5]. 

Second, let us explain the necessity of the regulator. 
Without it w e wo uld not have generated the first term on 
the RHS of (2.7), which eventually gives the desired re- 



sult. Instead, the other two terms on the RHS would have 
contained singularities for A„ + z — 0. These singularities 
need regularization. Note also th at th e regulator cannot 
be removed after the derivation of ( 2.7 ) since the first term 
on the RHS would then become ( U„ 7r6 2 (A„ + z )) Nf=2 - 
While this looks like the desired delta function we have to 
remember that for e = the fermionic measure is zero for 
A„ + z = 0, see (2.6). Consequently, at any finite volume 
this term is equal to zero for any z. Taking the thermo- 
dynamic limit does not make it nonzero. In conclusion, 
we only obtain the desired result by using a regulator and 
removing it at the very end of the calculation. This is in 
contrast to the derivation of the usual Banks-Casher rela- 
tion at fi = 0, where the factors of A n + m in the fermionic 
measure are always nonzero since in this case A„ is purely 
imaginary and m is real. 

Third, we comment on the possibility of choosing real 
instead of complex conjugate masses. With m\ — m 2 = m 
the RHS of ((O) becomes (3A 2 /tt 2 + 2H 2 )m 2 so that we 
cannot extract A 2 alone by taking derivatives with respect 
to mp^j However, as an alternative to ( 2.5 1 we can compute 

Id 2 3 
lim lim lim — — log Z(m\, m 2 ) — ~^A 2 

m-»OV r 4->oomi 1 2=m V4 Om\Om-2 7T Z 

(2.14) 

to extract A 2 . In the calculation of the other side of the 
equation we first compute d mi d m2 log Z(mi, m%; e) at fi- 
nite volume and then set mj = m 2 = m. After this step, 



11 There are typos in the spectral density in |41l eq. (38)] and 
H3 eq. (34)]. They must be divided by it. 

12 At lower density, where Ua(1) is explicitly broken by the 
axial anomaly, a similar problem occurs even with complex 
quark masses. In this case the coefficient Ci ns t mentioned in 
footnote [7] is nonzero so that derivatives with respect to z and 
z* extract a linear combination of A 2 and Cj ns t- 



z on the RHS of (2.7) is replaced by m and the measure 
is again positive definite. The rest of the calculation then 
goes through in the same way and yields ( |2.13 ). Never- 
theless, we have chosen to work with complex conjugate 
masses since in this case the measure is positive definite 
at all stages of the calculation, which is mathematically 
more rigorous. 

Finally, we discuss an alternative way to regularize 



(2.2) without diagonalizing D: 



Z(z, z*;j) = (det [(D + z)(rf + z*) + f] ) YM 
3 D + z 1 ' 



det 



(2.15) 



YM 



The regulator j has a clear physical meaning: It is a source 
for the diquark condensate in two-color QCD (and for the 
pionic condensate in QCD with nonzero isospin chemical 
potential, see section [3]). Elementary calculations lead to 



ddlog Z(z,z*-j) = 
1 

tr 



1 



(D+z){D^+z*)+f J (Dt + z*)(D + z)+j 2 
1 



tr 



x tr 



tr 



N f =2 



(Di + z*) 
(D + z) 



(D + z)(£»t + z *) + f 
1 



x ( tr 



(DUz*) 
(D+z) 



(£>t + z*)(D + z) +j 2 
1 



N f =2 



{D + z)(D^ + z*)+j 2 
1 

(Di + z*)(D + z)+j 2 



N,=2 



Nf—2 



(2.16) 



If the chiral limit z — > is taken at finite volume, the 
second and third line vanish because of chiral symme- 
try. Furthermore, it was shown in appendix A of [39] 
that at finite volume the first term on the RHS becomes 
(tt J2 n ^ 2 ( z + / ^"))a' / -2 m tne 3 ~* ® limit. However, tak- 
ing j — > before V4 — > 00 does not give the desired result, 
for the same reason that e — > before V4 — >■ 00 did not, 
see the second comment above. Unfortunately, we have 



been unable to rewrite (2.16) in terms of the eigenvalues 



of the Dirac operator D du e to i ts non-Hermitian nature 
and therefore cannot show (2.13) explicitly. This does not 



imply that the correctness of the latter should be do ubted . 
We just cannot provide an alternative derivation of (2.13) 
from (|2.16|). 



2.2 Consistency with the microscopic limit 

In [15] we proposed that the microscopic spectral density 
defined bJ^l 



(2.17) 



13 Because of the change in conventions compared to [T3] (see 
footnote 10 1 there is now no factor of I/V4 in front of pi. 
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is universal and can be computed in random matrix the- 
ory (RMT). This conjecture was consolidated in [TB] by ex- 
plicit construction of the RMT corresponding to two-color 
QCD at high density. The function p s (z) itself was com- 
puted in [17] for arbitrary even Nf. The explicit formula 
for p s (z) in the limit of maximum non-Hermiticity reveals 
that it asymptotically converges to an ./Vy -independent 
constant, 



lim p s (z) = - 

\z\— >ca TT 



(2.18) 



where we have to exclude purely real and purely imaginary 
values of z. This implies for large z 



SA 2pl 



1 

7T 



(2.19) 



Note, however, that the argument of p\ goes to zero since 
the limit V4 — > 00 is taken before the limit \z\ — > 00. Hence 



:T M0) = i 



3A 



(2.20) 



which exactly agrees with (2.13). Although this "deriva- 



tion" is by no means rigorous, it is reassuring to see that 
the microscopic limit is consistent with our main result. 



evaluated in the limit rnxmiV/^A — > 00, we find that the 



leading-order terms on the RHS of (2.231 cancel and the 
next-to-leadi ng-or der term gives ZViA z /it 2 , which is con- 
sistent with (2.21 ) . This consistency check is based on the 



fact that the e- and the p-regime have overlapping domains 
of validity 03111]. 

Let us now consider mj = m 2 = m. From (2.4) we 
then obtain 



lim lim — ^ 

ra-fOVi^oo K4 d(m 2 ) 



log Z{m, m) 



3 A 2 

TT<& 



2H 2 



while from (2.22) with mi = to 2 = m we obtain 



d 



d(m 2 ) 



log Z(m, m) 



-X 2 + m" 2 



(2.24) 



(2.25) 



N f =2 



To evaluate the RHS we use jT5j eq. (4.32)], evaluated in 
the limit rr^ViA 2 — > 00, and obtain W4A 2 /ir 2 . This is 
not consistent with ( |2.24[ ) since the H2 term is missing. 
The reason for this discrepancy is that in [15 the high- 
energy term H2 tr(Af 2 ) was not considered. Including this 
term modifies some of t he sp ectral sum rules and makes 
(2.25) consistent with (|2.24|). This issue is discussed in 



more detail in appendixTB 



2.3 Consistency with spectral sum rules 

We now consider real quark masses. If we are not inter- 
ested in deriving a Banks-Casher-type relation there is no 
reason to insist on a positive measure, and we can take 
mi m 2 - We also omit the regulator e. 

From (2.4) we obtain with z — > mi and z* — ► m 2 



i,m 2 -v0 V 4 ->oo V4 omidm2 ft 



(2.21) 



On the other hand, starting from the first line of (2.2) 



and using the fact that the Dirac eigenvalues come in pairs 
±A n due to chiral symmetry, we have 

Z(mi,m 2 ) = (l£(-\l + mi)(-\ 2 n + m 2 2 )) , (2.22) 

\ n I YM 

where the prime on the product (and on the sums below) 
means ReA„ > 0. This gives 



d 2 log Z 
dmidm 2 



^— y' 2mi -i ^' 2m 2 

2 1 m 22j 



ni 

2mi 



-\ 2 m + m 2 ^ -X 2 



N f 



2m 2 



-A 2 + m\ 1 Nf=2 



(2.23) 



To compute the RHS we can make use of the massive 
spectral sum rules derived in |15j in the e-regime of chiral 
perturbation theory. From eqs. (4.32) and (4.33) of that 
reference (with A — > iX due to the change in conventions), 



2.4 General Nf 

So far we have worked with two flavors, and we now show 
that our main result is unchanged for general even Nf. 
We have nothing to say about odd Nf since the measure 
is not positive definite in this case. The mass matrix is now 
M = diag(z, . . . ,z,z* , . . . ,z*) with Nf/2 entries each of 
z and z* . T he ge neralization of (2.3) to even Nf > 4 is 
given by [15j^ " 



U 15 



5£ 



3A 2 

47T 2 A^SZS R 



max Re\A 2 tT(MS R M T Sl) 



A* 2 tr(MA7 L A/ T A7]j)| , (2.26) 



14 See the discussion before (2.3). In |15| we assumed Mlr = 
M and M RL — M\ while we now have Mlr — Mrl = M. 
This implies that additional terms in the chiral Lagrangian 
are allowed by symmetries, e.g., A 2 tr(M ErM* However, 
HDET shows that the coefficients of these terms vanish and 



that only the terms in 1 2.26 1 remain. See also |45II46| for an- 
other example where terms allowed by symmetries vanish in 
HDET. Note also that for the values of Er and El that max- 
imize ( |2.26[ ) (see below) the term in ■ ■ ■ ■ } is already real. 
15 Similarly to footnotes [7] and 12 for Nf = 2, the axial 
anomaly (or instantons) gives rise to another mass term at 
0(rn 2 ) in the free energy for Nf > 4, which is not displayed 
in (2.26 1 because it vanishes at sufficiently high density. But it 



is nonzero at lower density. This mass term is obtained from 
the instanton-induced interaction with Nf left-handed and Nf 
right-handed fermion legs, by replacing two sets of left- and 
right-handed fermions by M 2 and the remaining 2(Nf — 2) 
fermions by diquark condensates, in a way analogous to |19l 

m. 
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where A £ U(l), E L <= SU(A^ / ) L /Sp(7V / ) L and (L <H> R). 
An explicit parameterization of E L /r is given in |15j , but 
it is unimportant for our present purpose. Note that the 
minimum has become a maximum together with an overall 
minus sign. We now use a Cauchy-Schwarz inequality for 
matrices, 



\ti{X^Y)\ 2 < tr(X t X)tr(Y t F) 



(2.27) 



which becomes an equality if and only if X and Y are 
multiples of each other [37] • Because of E l El = Er^Jr = 
1 this leads to 

\/r T y?t 



| tr(ME R M T El)\ < tr(AftM) - JV/zz* , (2.28) to obtain Kl2h 



3 QCD at nonzero isospin density 

Let us now consider QCD with two flavors and any number 
of colors N c > 2 at large isospin chemical potential \ii = 
2\x. In this case we have 

Z(m ll m 2 ) = (det(D(p) + m x ) det(D(-(i) + m 2 )) YM 

= (det(D + m 1 )dct(^ t +m 2 )) YM , (3.1) 

where the last line follows from D(fiy = —D(—fi) and the 
fact that the eig envalues of D come in pairs ±A. This is 
the s ame as (2.2). Setting M — diag(z, z*) we again s tart 
from (2.6) and go through the same steps as in section 2.1 



where for our choice of M the maximum is obtained^ for 
E L = Eii = H with 



I = 







L N f /2 





(2.29) 



Together with A = ±1 we thus have 
58 



2tt 2 



(2.30) 



Including the high-energy term H 2 ti:(M 2 ) we obtain 



1 N f 

- log Z(z, z*) = — 



3A 2 



zz*+H 2 (z 2 + z* 2 ) 



■0(\z\ 3 ). 
(2.31) 



Thus the only diff erence to (2.4) is a factor of AT//2, and 
the analog of (2.5) is 



lim lim — ddlos Z(z, z* , 
z-s-o v 4 -s-oo V4, 2tt 2 



m *A 2 . 



(2.32) 
The analog of ([2^ is 

Z(z, z*; e) = ( J] [(A„ + z)(X* n + z*) + e 2 } Nf/2 ) yM , 

n 

(2.33) 



and going through the same steps as in section |2.1| we 
obtain 



lim lim lim — dd log Z(z, z* : e) 

e->0 2->0 y 4 ->-oo V A 



N 



(0) 

(2.34) 



with the one-point function in the presence of Nf massless 
flavors. Hence the Banks-Casher-type relation for Nf > 4 
is the same as for two flavors. The arguments of sec- 



tions 2.2 to 2.3 also generalize straightforwardly to general 
even Nf. 



16 Note that the maximum could also be obtained for El = 
1 N 



Er = ± ^ J, but such E are not elements of the coset space 
SU(JV»/Sp(JV» and therefore excluded. 



In QCD at sufficiently large fij, BCS pairing of the 
form (dr/su) is formed near the Fermi surface according to 
the BCS mechanism 23 , 24J . The chiral Lagrangian in this 
regime can be constructed similarly to 15J, but the coset 
space is now 11(1)^4 x U(l)/ 3 [IE], where the latter is the 
U(l) symmetry with respect to the third isospin genera- 
tor. Repeating th e an alysis in 15j with the modification 
discussed before (2.3), the mass term can be written as 
(4.13)] 



eq. 



c iso Z\ 2 {(v4 2 + ^* 2 )detM}, (3.2) 



where A e U(l) and c- lso is a positive coefficient, see (3.10) 
below. Thus 



5£ = — Ci SO A 2 max 



{(A 2 + A* 2 ) det M } = -2c iso A 2 



(3.3) 



where the last equation is obtained with A = ±1. 

To compute the coefficient c- lso we use essentially the 
same arguments as for the color-flavor-locked (CFL) phase 
[35] and for dense two-color QCD [15], the only new fea- 
ture being that the sign of the chemical potential is op- 
posite for the two flavors. A simple way to adapt QCD 
with isospin chemical potential to the method of [551IT5] 
is to change the basis of the quark fields from (u, d) to 
(u,d c ) with d c = Cd T . Accordingly, the Lagrangian is 
then rewritten as 

C = u[D{ji) + TOi]tt + d[D(— ji) + m 2 ]d 

= u[D(fj) + mi]u + d c [Z?af(m) + m 2 ]d c , (3.4) 

where -Daf is the Dirac operator in the anti-fundamental 
representation of SU(A r c ). Adapting the mass term in 
HDET, see, e.g., [T5] eq. (A.l)], to the present case we 
obtairO 



-HDET 
■'mass 



8/i 4 

x {f A ) ac {T A ) hd M lk M 3l + (L+*R), (3.5) 



17 The notation in ( |3.5[ ) is a bit sloppy since C is a 4 x 4 ma- 
trix, while the right- and left-handed spinors are two-spinors, 
which here are considered to be four-spinors with two compo- 
nents set to zero. Choosing C 

(4)" - ■ '' 



17472 we have {dl) a L Cu\ 



la 2 u L = (4°t£)* and {d T c Y R Cu d R -> -(df)%a 2 u d R 
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where ipi = u and "02 = d c , g is the coupling con- 
stant, a, b, c, d are color indices, k, I are flavor indices, 
and the f A = -(T A )* are the generators of SU(iV c ) in 
the anti-fundamental representation with normalization 
tr(T A T B ) = 5 AB J2. To leading order in g 2 , the expec- 
tation value of (3.5| factorizes into the product of the 



expectation values of the diquark condensates, and the 
latter can be computed in a weak-coupling calculation 
[49] . At high isospin density, BCS-type diquark pairing 
occurs in the color-symmetric, flavor-antisymmetric, and 
pseudoscalar channel, i.e., the quantum numbers of the 
condensate are the same as for pions [23l24| . We therefore 
have 



<(</>DiC(rf) 



3/N f =2 



S ab eijr L A L and (Lai?), 



(3.6) 



where r L / R G H>o an d A L / R e U(l) are magnitude and 
phase of the pionic condensate. For the magnitude we can 
use the result 50, cq. (3.42)], which in the present context 
reads 



(d T c C lb u) Nt=2 | = 2 



/ 6iV c n 2 A 
N 2 - 1 ~^g~ 



(d T C Cu) N}=2 | = 



This implies 



r L = Tr 



I 6N C n 2 A 
N 2 ~ 1 ~^g~ 



(3.7) 



(3.8) 



-(N 2 — 1) we thus obtain with 



Taking the expectation value of (3.5) and noting that 

gab ficd ^rpA^ac ^pA\bd 

U = A\A R 



5£ 



3N, 



c A 2 max{(L/ + C/*)detM} 



37V, 



C A 2 ~~* 



zz 
(3.9) 



with U = 1 in the last step. Comparing with (3.3) we find 

(3.10) 



3N n 



ISO ~ 4tt 2 ' 

Including the high-energy term H2 tr(M 2 ) we thus have 

y log Z(z, z*) = ^A 2 zz* + H 2 (z 2 + z* 2 ) + 0(\z\ 3 ) . 

(3.11) 



From (2.12) and (3.11) we obtain the Banks-Casher-type 
relation 



(3.12) 



For N c — 2 ( 3.12 ) reduces to ( 2.13 ) as it should (recall that 
there is no difference between isospin chemical potential 
and baryon chemical potential for N c = 2 51j). 



4 Adjoint QCD 

We now derive the Banks-Casher-type relation for QCD 
with fermions in the adjoint representation. To obtain the 
Dirac spectrum in adjoint QCD one might naively want to 
consider the Nf — 1 partition function (det(D + z)), with 
D in the adjoint representation of SU(A r c ), but for our 
purposes this is not a sensible choice because the positivity 
of the measure is spoiled for complex z[^]For this reason 
we consider an even number Nt of flavors for which we can 
insert complex conjugate masses, i.e., the mass matrix is 
M = diag(z, . . . , z, z* , . . . , z*) as in section [2~4] and 



YM 

(4.1) 

where the last equality follows from the anti-unitary sym- 
metry {Cj$K, D] = of the adjoint Dirac operator [5], 
which implies D* = (fr y 5 )- 1 D (Cm). On the QCD side 



Z(z, z*) = Uet N ' /2 {D + z) det N f /2 {D + z*)) 
, det N ^ 2 (D + z)det N ^ 2 (D^ +z*)) 



we are then led to (2.34) as described before for two-color 
QCD. 

We now turn to the HDET side. The construction of 
the chiral Lagrangian proceeds in close analogy to (15 , 
the difference being that the coset space is now U(l)s x 
JJ(1) A x [SV(N f ) L /SO(N f ) L ] x [S\J(N f ) R /SO(N f ) R ] [TS]. 
Repeating the analysis in (15 with the modification dis- 
cussed before we obtain HfiP 1 



= -c adi A 2 {A 2 tv(MS R M T sl) 

+ A* 2 tr(M£ L M T £ R )} , (4.2) 

where A e U(l), S L € SU(iV / ) L /SOnV £ ) L and (L a R), 
and c a( jj is a positive coefficient, see (4.8) below. Thus 



5£ = -c adj A 2 max Re \ A 2 tv(MS R M T s[) 

+ A* 2 tr(MS L M T r]j)| 
= -c adj Z\ 2 -2N f zz* , (4.3) 



where we have employed the argument leading to ( 2.28 ) 

The coefficient c a( jj can again be computed From 
HDET. Instead of (3.5) we now have 



C 



HDET 



^((4)?C(^)i)((^)fcC(^)f 



jr 

8^' 

x (^) oc (t A ) bd M ife M, 7 + (L a R) , 



(4.4) 



18 For real mass mass z — m the measure is actually pos- 
itive in Nf — 1 adjoint QCD so that in principle one could 
derive a Banks-Casher-type relation. However, in this case the 
terms involving A 2 and H2 have the same mass dependence so 
that A 2 alone cannot be extracted by taking derivatives with 
respect to m. 

19 As in footnote 



14 



additional terms are allowed by symme- 
tries, but their coefficients are zero in HDET. 
20 The maximum is now obtained for A = ±1 and El = 



1 

Nj/2 



N f /2 
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where the t A (A = 1, . . . , JV| — 1) are the generators of 
SU(iV c ) in the adjoint representation, i.e., (t a ) bc = —if abc 
with the structure constan ts f abc . To leading order in g 2 
the expectation value of (4.4) again factorizes into the 
product of the expectation values of the diquark conden- 
sates. The condensation again occurs in the pseudoscalar 
channel, but another difference to [T5] is that the diquark 
condensate is now symmetric in the color and flavor in- 
dices (and antisymmetric in the spin indices as before), 
i.e., 

(^'E)tC^ L ) b j ) Nf =5 ab 5 ij r L A Li and [L R) , 

(4.5) 

where r L / R € R> and Ar^/mi £ U(l) are the magni- 
tudes and phases of the condensates. From weak-coupling 
calculations at high density we fincj^] 



tl = r R 



3 [i 2 A 



N c ng 



(4.6) 



While the magnitude of the condensates is fixed regardless 
of the mass term, the phases of the condensates must be 
such that the ground state energy is minimized. Taking the 
expectation value of (El and using S ab 6 cd (t A ) ac (t A ) bd = 
-N C (N 2 - 1) we thus obtain 



se = - 



3(iV c 2 - 1) 
8tt 2 



A 2 



x max 

A 



ax Re { ti(A L MA R M T ) + (R o L)\ 



3(iV c 2 - 1) 
8tt 2 



A 2 ■ 2N f zz* 



(4.7) 



where Al — diag(Aki, . . . , Alat,) and (L O R), and in 
the last s tep we have once again used the argument leading 
to (2.28). Note that the maximum is now obtained for 
A R = M~ x A h M* pi Comparing with dOl we find 



Cadj 



3(jV c 2 - 1) 
8tt 2 



(4.8) 



Including again the high-energy term H 2 tr(M 2 ) we have 
1 



Va 



log Z(z,z*) 



3(iV c 2 - 1) 
4tt 2 



N f A 2 zz* + H 2 (z 2 + z* 2 ) + C(|z| 3 ) . (4.9) 



Comparing with (2.34) we arrive at the Banks-Casher- 
type relation for adjoint QCD, 



A 2 = 



2 ^ 3 



3(iV 2 -l)' 



(4.10) 



21 This can be obtained by replacing the group theoretical 
factor in the gap equation in |351|50] by = 

-N c {l) be . 

22 This assumes that M is diagonal since for general M the 
combination M~ AlM* would not be diagonal. 



Let u s ad d an important remark here. In deriving the 
RHS of (4.7) it was essential to include the phases of the 



condensates. The minimum of 8£ is obtained only if some 
of the phases are nontrivial, i.e., 5£ would not have been 
minimized if all phases had been set to unity. In contrast, 
in the other two cases (two-color QCD and isospin) the 
minimum of 8£ is obtained for trivial phases, but this 
should better be regarded as accidental. 



5 Conclusion and outlook 



In this paper, we have derived a new Banks-Casher-type 
relation at high density which relates the complex Dirac 
eigenvalue density at the origin to the BCS gap squared. 
We have shown that the proportionality constant in this 
relation is independent of the number of flavors and is 
fixed by weak-coupling calculations. (The latter point 
should be contrasted with the usual Banks-Cashcr relation 
in the QCD vacuum, where the proportionality constant is 
fixed by the definition of the chiral condensate.) We have 
also checked the consistency of our relation with previous 
results for the spectral density of Dirac eigenvalues in the 
microscopic limit [T5"|I17|. 

As mentioned in the introduction, the positivity of the 
measure is a necessary condition to derive the Banks- 
Casher-type relation. This is why its generalization to 
the CFL color superconducting state [S5] in QCD at high 
baryon density is not feasible, though the Dirac eigenvalue 
spectrum itself is related to the BCS gap [H]. Still, it is 
likely that the BCS gap of the CFL phase is characterized 
by some singular behavior of the Dirac spectrum near the 
origin, like the scenario in [27U53) . This would be an in- 
teresting question to explore. 

Our relation enables us to compute the magnitude 
of the BCS gap numerically through the low-lying Dirac 
eigenvalue spectrum on the lattice. More detailed informa- 
tion on the Dirac spectrum is provided by the microscopic 
spectral density (and higher-order spectral correlations) 
of the Dirac eigenvalues at high density. These quanti- 
ties have been computed for two-color QCD [T7], and it 
would be interesting to compute them for the isospin and 
adjoint cases as well. Combined with the measurement 
of the diquark (or pionic) condensate that can be ob- 
tained from the Dirac singular value spectrum |18j . one 
should be able to perform a detailed quantitative analysis 
at high density numerically. We hope that such calcula- 
tions from first principles will eventually lead to solid ev- 
idence that the BCS mechanism, and moreover the BEC- 
BCS crossover phenomena, ubiquitously occur not only in 
condensed matter systems but also in high-energy physics 
(QCD). 
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A Computation of the high-energy constant 

Hi 



The term H 2 tr(M 2 ) in ( pH| and ( |2.31| ) was originally in- 
troduced in chiral perturbation theory as a counterterm 
needed in the renormalization of one- loop graphs |54[I37] . 
For the same reason it is necessary in the chiral pertur- 
bation theory of two-color QCD at high density. At suffi- 
ciently high density where g <C 1, the coefficient H 2 can 
be determined by computing the ground-state energy den- 
sity microscopically. In this appendix, we compute H 2 to 
leading order in g (i.e., in the free theory). For this pur- 
pose, we use the massless quark propagator with chemical 
potential /up*] 



where we have regularized the divergent integral by in- 
troducing a momentum cutoff ^4tjv- Note that this term 
depends on the scheme used to define the finite part of 
the diagrams. For instance, dimensional regularization or 
Pauli-Villars regularization would lead to different values. 
If we had used dimensional regularization the divergent 
integral in (A.4) would simply be zero. 



Putting the terms together and comparing with (2.31 ), 
we obtain 

H^-^^+A^), (A.5) 



4tt 2 



where the second term depends on the scheme we employ. 
This result is also valid for the other cases (QCD with 
isospin chemical potential and adjoint QCD). 



1 



(Po + M)7o - P • 7 



f + mo [po + Ai + iesgn(p )] 2 P 2 

= {Po + m)7o -P I 
Po+ fJ.+p — ie 

0(p-n) 6{H-P) 



Po - (p- M) + i£ Po - (P - A 1 ) - 



(A.l) 



where sgn(x) denotes the sign of x and p = \p\. A one-loop 
calculation gives the free energy due to the (degenerate) 
quark masses, 



Z(m) 

W) 



(-1) 



V A m< 



d 4 p 



Tr 



1 



f + ^lo 



, (A.2) 



where the factor (—1) is due to a fermion loop and "Tr" 
is the trace over spinor, color, and flavor indices. When 
the integration contour is closed in the upper half-plane, 
the integrand has two poles at pq = ±p — \i. From ( A.l ), 
the pole at po = p — \i gives rise to the /x-dependent con- 
tribution, while po = —p — fi gives the /i-independent vac- 
uum contribution. Performing the integration over p, the 
/^-dependent part reads 



-d Tep N f V 4 



d 3 p 
(2tt) 3 

2 2 
\x m 

4tt 2 ' 



tr 



d 4[(p 
dpo (Po 



H) 2 -P 2 
p + Id) 2 



Po=p-n 



(A.3) 



B Modified spectral sum rules 

The purpose of this appendix is to discuss modifications 
of the spectral sum rules due to the high-energy constant 
H% that was not taken into account in |15j . For simplicity 
we focus on two-color QCD at high density with Nf = 2 
and real masses. In the e-regime at high density we have 
from Q and [B] 



Z(mi,m 2 ) = Io(^ViA 2 mim2) exp \V i H 2 {m 2 l 



"■2)1 ' 

(B.l) 



where I a is a modified Bessel function and the second fac- 
tor is missing in [TS]. This is the most generic form of 
the partition function at leading order of the e-expansion, 
consistent with all symmetries of the system. 

We now investigate how the inclusion of the high- 
energy term modifies the spectral sum rules derived in 
|15j . Taylor expansion of (B.l) gives 



Z(mi,m 2 ) = 1 + V 4 H 2 (ml + ml) + ot(y 4 A 2 ) 2 m\m\ 

+ hy 4 H 2 ) 2 (m\ + ml) 2 + ^{V A H 2 f{m\ + mj) 3 

+ a(y 4 A 2 ) 2 (y 4 H 2 )m 2 m\{m\ + m 2 2 ) + 0(m 8 ) (B.2) 

with a — 9/(47r 4 ). In terms of the complex Dirac eigen- 
values A„ the partition function can also be expressed as 



where "tr" is now the trace over color and flavor indices 
and d rep is the dimension of the fermion representation. 
On the other hand, the vacuum contribution is given by 



£ M =o(m) 



-V 4 m 2 
2 



d 4 p 



Tr4 



-d Icp N f V 4 ATy 



rn 



Air 2 



(A.4) 



23 Note that the contour in the usual Feynman prescription is 
modified from fj, = to fj, / [55l fig. 4]. For fj, = 0, ( |A.1| ) re- 
duces to the conventional is prescription because po sgn(po) > 
0. 



Z{m u m 2 ) = l\{(l- r f] (l 



N f =2 



2 \2 



\k<i 



k /y l 



1 2 2 

+ m 1 m 2 



\ n A "/ / \k<£< 



I 1 



i\m 2 2 (m\ + m 2 ) ( E't-T E' lAa > + ( m ") 
\ n A « k<e X k\ I 



(B.3) 
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where the prime on the sums (and the products below) 
again means Re A„ > 0. Starting in the second line we 
ha ve om itted the subscript Nf = 2 for brevity. Matching 
of (B.2) and (B.3| yields massles s sp ectral sum rules that 
incorporate the effect of H 2 , e.g.j^j 




V 4 H 2 



a(V 4 A 2 ) 2 + [V 4 H 2 f 



( \ - \2\ 2 

\k<Km A k A e A m t 

yV 1 yV 1 

\ n A « k <e X k X i , 



—x{ViH 2 ) 




(B.4) 
(B.5) 

(B.6) 
(B.7) 



{V 4 H 2 f 



a{V 4 A 2 ) 2 (V 4 H 2 ), 

(b.s; 



Thus the spectral sum rules are drastically modified after 
the inclusion of the H 2 term. However, it is intriguing that 
we can nonetheless find special combinati ons o f sum rules 
that do not depend on H 2 . For instance, (B.6) — 2 • (B.5 1 
yields 




a(V 4 A 



2\2 



(B.9) 



which agrees with the result of [TS]. It seems likely that 
there exist infinitely many such spectral sum rules that 
are not affected by the H 2 term. If true, the microscopic 
spectral correlation functions of RMT that ignore the H 2 
term still describe at least some portion of the universal 
spectral correlations of the small Dirac eigenvalues. In the 
rest of this section, we show that the correlation functions 
of Dirac eigenvalues of order \/^JV 4 A 2 are still described 
by RMT and unaffected by the H 2 term. 

To make sense of the modified sum rules, let us assume 
that the Dirac spectrum {A„} is a superposition of two 
statistically independent spectra {a„} and {A n } such that 
{A„} = {a n } U {A„}. Then the expectation value in the 



chiral limit (• • •' 



is understood as 



(0) Nf - 



' „2T1' A2\ N f 



[da][dA]0(U'AU' t A 
[da][dA](U' k alU.iAj) N ' 



(B.10) 



It is essential that the measures for {a n } and {A n } are 
factorized, i.e., that the two spectra are distributed inde- 
pendently. 

Moreover, we assume that the eigenvalues a n are of 
order 1/ ' \JV 4 A 2 and satisfy the spectral sum rules without 



24 Similarly, we could derive modified massive spectral sum 
rules following the approach of [151 sec. 4.3]. 



the H 2 term, i.e., 



IT 



1 - 



1 - 



;V 4 A 2 miiri2 



(B.ll) 

whereas the A n are solely described by the H 2 term, i.e., 



ir 



1 



A^, 



1 




exp \y 4 H 2 (m\ + ml)] , 
(B.12) 



where we have again omitted the subscript Nf — 2 for 
brevity. Under these assumptions, the mass dependence 
of the partition function is reproduced correctly since 




A 2 



= I -^V 4 A 2 m 1 m 2 exp \V 4 H 2 (m\ 



(B.13) 



It follows from this that all spectral sum rules for {A„} can 
be deduced from those for {a n } and {A ra }. For instance, 




2 • • {~V 4 H 2 ) + (v 4 H 2 y 



a(V 4 A 2 ) 2 + (V 4 H 2 



(BAA) 



which agrees with (B.6). The other sum rules can be re- 
covered similarly. 

In our present argument, the spectral correlation func- 
tions obtained from RMT are exact results for {a n }, but 
not for {A„}. At this stage we know nothi ng ab out {A n } 
except the mass dependence postulated in (B.12). We now 



propose that {^4„} is the spectrum of the Dirac operator 
in the free theory, £>f rcc = 7iy<9„ + njo- This may be un- 
derstood as follows. From the result obtained in appendix 
[A] we have 



log 



Z(m) 

W) 



= V 4 [H 2 (m 2 +m 2 ) +C(m 4 )] , (B.15) 
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where there is no C(m 3 ) term on the RHS because of chiral 
symmetry. In the e-regime at high density, where U 4 — > 00 
with rn ~ 1 / ' \JV^A 2 , the higher-order terms in m drop 



out, and we obtain log[Z (m) / Z (0)] 
an exact result, i.e., 



V 4 H 2 (m( 



Z(m) 



= det' 1 - 



n 2 



1 - 



D 2 



= cxp [V A H 2 {m\ + ml) 



(B.16) 



There fore the eigenvalues of the free Dirac operator satisfy 
(B.12), and the {A n } are scattered over the entire range 
from 1/L to \JH 2 , where L is the linear extent of the 
box. On the other hand, the typical scale of the {a n } is 
\j\JV 4 A 2 1 which is much smaller than 1/L. This means 
that the domains of {A n } and {a n } do not overlap. The 
existence of the H 2 term does not change the spectrum on 
the scale 1 / \/V 4 A 2 at all, and our previous results [TolTT] 
based on RMT remain exact as far as the eigenvalues on 
this scale are concerned. 

Finally, let us recall the e-regime at /i = [5] for com- 
parison. There we have no H 2 tr(M 2 ) term because it is 
subleading in the e-expansion. What happens if we add 
this term to the chiral Lagrangian by hand? Then the 
leading spectral sum rule will be altered to 



E 



> 1 

A2 



diViE) 2 + C 2 V 4 H 2 (B.17) 



with numerical constants C\ and C 2 . The second term may 
be interpreted as the contribution from the perturbative 
Dirac spectrum since 




free 



W „ V4A 3 

dX — 



uv ■ 



(B.18) 



However, in the limit V4 — > 00 with A\jy fixed, the first 
term in (B.17 1 dominates the second term, and a con- 
ventional Leutwyler-Smilga sum rule is recovered^] The 
contribution of small eigenvalues dominates the UV part. 

On the other hand, the contribution of small eigenval- 
ues in dense QCD is not large enough to overwhelm the 
UV part. For example, in (B.4) the small eigenvalues con- 



tribute nothing, while the U V s pectr um yields V 4 H 2 in 
the same manner as in (B.18). In 



interacting case pi(0) is nonzero due to the spontaneous 
breaking of chiral symmetry. To extend this comparison to 
nonzero chemical potential we now compute the spectral 
density of the free Dirac operator at /i ^ 0. To the best of 
our knowledge this is a new result. 

The free Dirac operator in Euclidean space (where we 
assume even dimension d) reads 



D (v) = (7>A + M7o) ® ld rcp with v = 1, 



(C.l) 



where d rep denotes the dimension of the color representa- 
tion in which the quarks transform if coupled to the gauge 
field. The spectral density of D{p) is defined asQ 

Pl (z = x + iy) = ^-txS 2 (z-D(fi)), (C.2) 
Vd 

where Vd = L d is the volume of a (i-dimensional torus with 
linear extent L. Our main result is 

Pl(*) = drep^^^0(M-M) 

X(x 2 + y 2 )[(n 2 -x 2 )(n 2 + y 2 )}^ , (C.3) 

where C D = 2n D/2 / T(D /2) is the area of a D-dimen- 
sional unit sphere. Note that (C.3) is a continuum result 



and does not reproduce the spectrum of the free lattice 
Dirac operator. It is easily seen that our result has the 
correct mass dimension (d — 2) and that pi(0) = in the 
free case. Note also that for \z\ <C the effect of fi formally 
drops out and pi{z) cx \z \ 2 . S napshots of p\{z) are given 



for d — 2 and 4 in figure C.l (where all parameters were 
made dimensionless by an arbitrary refe rence scale). 
Let us also check the /i — > limit of (C.3). Using 



lira 



(M 2 



,d-2 



■0fr-\ X \) = ^ r «*7* )/2) 8( x) 



r(d/2) 



(C.4) 



we obtain 



lim pi(z) 
fj,— >o 



(2Tr) d / 2 r{d/2) 



{yl^Six), (C.5) 



of the same order in the thcrmoi 
cannot drop the UV piece alone. 



RGj) the two pieces are wn j cn j s tne correc t free spectral density at /z = 0, see 



ynamic limit, and one 



for d rcp = N r = 3 and d — 4. 



C Dirac spectrum in the free limit 

At zero chemical potential there is a marked difference be- 
tween the free and the interacting Dirac spectrum. In the 
free case the spectral density behaves like A 3 , while in the 



We now present the derivation of (C.3). Using the re- 
lation ^(l/z) = ttS 2 (z) we can obtain the spectral density 
from the resolvent of -D(/i), which is defined by 



G(z) 



V d z 
d 2 d/ 



1 



2 / d d p 



25 It is essential that the cutoff is fixed in the thermodynamic 
limit. In actual lattice QCD simulations, the cutoff (the inverse 
lattice spacing) and the volume are finite, implying that the 
Leutwyler-Smilga sum rules will not be accurately satisfied. 



(27r) d z 2 + (p + iii) 2 +p 2 



(C.6) 



26 Whether or not the factor of 1/Vd is included in the defi- 
nition is a matter of convention. Here we are consistent with 



12.8) 
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Fig. C.l. jOi(a) for d = 2 (top) and d — 4 (bottom) at fi = 1. 



where we have used the fact that the eigenvalues come in 
pairs ±A, p denotes a (d — l)-dimensional spatial vector, 
and the factor of 2 d / 2 is due to the 7-matrices having 
dimension 2 d / 2 . This gives 



Pl (z) = -dG(z) 

7T 

2 d/2 f d d p _^ 



0- 



rcp 7T J {2n) d "z 2 -w 2 

with w = i\J (pq + ip) 2 + p 2 
2 d/2 f d d p 



' • 2 J (2tt) 
2 d ' 2 C d ^ 



j[6 2 (z-w) + 6 2 (z + w)] 



rcp 2 (2tt)<* 

JdpoJ dRR d - 2 [5 2 {z-w) + S 2 {z + w)], 

(C7) 



where in the last line we changed variables to spherical 
coordinates with R = \p\. Writing w = wr + iwj we have 



wj — Wji = p^ — p + R and wrWi = —/J^po- Therefore 

1 
p 



R = Wj - w" R , — Pq + pf 
and thus 



2 (p 2 -w 2 R )(p 2 +w 2 ) >0 



(C.8) 



(C.9) 



The Jacobian of the transformation from (po,R) to 
{wr,wi) is given by 



d(p , R) 



Collecting all pieces we obtain 



pR 



(CIO) 



■2 d/2 gd-i 
2 (27r) d 



x / dw/i? d - 2 ^-[<5 2 (z-u;) + 5 2 (z + W )] 



2^^-! , 3 |z| 2 

v (M_ll) 



with R = - v/(/i 2 - x 2 )(p 2 + y 2 ) 
A 4 



(C.ll) 



which is ( C.3 ). 



D Replica trick 

The purpose of this section is to present an alternative 
way to derive the new Banks-Casher-type relation. We 
use the so-called (fermionic) replica trick. While this trick 
has frequently been used in theoretical physics over the 
years it can hardly be called mathematically rigorous due 
to a notorious ambiguity in continuing from integer n to 
real n, where n is the number of replicas. We refer the 
reader to [571H21I55] for detailed discussions of the pros 
and cons of the replica trick. In this section we do not 
discuss the validity of the replica trick but simply use it 
as an easy cross-check of our main result (implying that 
the replica trick actually works in the present case). To 
add some variety we perform the calculation for adjoint 
QCD rather than for two-color QCD. Our calculation here 
can easily be adapted to two-color QCD and QCD with 
isospin chemical potential as well. 

Since the Dirac operator at nonzero p is non-Hcrmi- 
tian, we must insert pairs of complex- conjugate replica fla- 
vors [59,31,60 . Denoting the number of dynamical flavors 
by Nf and the number of pairs of replica flavors by n, we 
have 



Pi (z) = 



1 

v 4 



lim — dd\ogZ N+2n (z,z*) (D.l) 



with 
Z 



Nf+2 



n (z,z*) = (det N f (D + m)|dct(£> + z)| 2 ™\ 



YM 



det JV ' (D + m)det"(D + z)det n (L> + z* 



YM 



(D.2) 
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where we used the property C75.DC75 = D* of the Dirac 
operator in the adjoint representation. To evaluate (D.2| 



at large \i we recall the shift (4.3) of the ground state 
energy density!^] 



se 



m 2 c 



8tt 2 



^A 2 max Re { A 2 tr {ME r M t E{) 
+ A* 2 tx{ME L M T E R )\ , (D.3) 



where the difference is that now Sl/r G S\J(Nj + 2n)/ 
SO(N f + 2n) and 



M = diag(m, . . . , m, z, 



, z*) with m real . 



N f 



(D.4) 



The maximum in (D.3) is attained with A = ±1 and 



E L = E r = ± 




(D.5) 



and thus 5£ 



3(N 2 



An 2 



-A 2 (N f m 2 + 2nzz*). Substitut- 



ing this and logZ = —V^SS into ( D.l[ ) we find 

] A 2 



Pl W ~ 27T3 



(D.6) 



This result agrees with ( |4.10 ) . 

By evaluating 58 to higher orders in M we could ob- 
tain a high-density analog of the Smilga-Stern relation 
[61] . However, this is beyond the scope of this paper and 
deferred to future work. 
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